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Abstract 

The thermodynamic singularities of QCD in the plane of complex baryo-chemical potential ^ 
are studied. Predictions are made using scaling and universality arguments in the vicinity of the 
massless quark limit. The results are illustrated by a calculation of complex ^ singularities in a 
random matrix model at finite temperature. Implications for lattice QCD simulations aimed at 
locating the QCD critical point are discussed. 
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I. INTRODUCTION 



The search for the QCD critical point has attracted considerable theoretical and exper- 
imental attention recently. The existence of such a point - an ending point of the first 
order chiral transition in QCD - was suggested a long time ago jl[ 01, and the properties 
were studied using universality arguments and model calculations more recently 0, 0] (see 
Ref. 0] for review). The experimental search for the critical point using heavy ion collisions 
has been proposed in 0]. It is apparent that theoretical knowledge of the location of the 
critical point on the phase diagram is important for the success of the experimental search. 

First principle lattice QCD calculations aimed at determining the location of the critical 
point on the T, /i (temperature, baryon-chemical potential) diagram have been attempted 
recently using several different techniques ^ ^ ^, J^] (see Ref. Oil for review). The major 
obstacle for direct Monte Carlo simulation is the well-known lack of positivity of the measure 
of the path integral defining QCD partition function at nonzero baryo-chemical potential /i 
— the sign problem. One of the methods to deal with this problem is to Taylor expand the 
QCD pressure in powers of /i around /i = 0, i.e., around the point at which direct Monte 
Carlo simulations are not hindered by the sign problem [sl. [lol|. 

The success of such an approach in determining the location {Te-, ^ie) of the critical ending 
point crucially depends on the convergence radius of the Taylor expansion around /i = 
HHEl. The convergence radius, in turn, is a function of the position of the singularities 



in the complex plane. Little is known about the location of these singularities to date. 
The purpose of this paper is to expand our knowledge of the location of these complex plane 
singularities. 

We shall be able to determine the position of the singularities in the regime where the 
quark masses m are sufficiently small. To summarize, the strongest rigorous consequence of 
our analysis is that the convergence radius, /i^, achieves its minimum value at a temperature 
slightly above Tc (by 0{m}/^^^'>)). This value scales as 

mm/ifi(r) ~m^/(2/35)^ 

and vanishes in the chiral limit (m 0). The value l/(/35) ~ 0.54 is determined by the 
critical exponents of the 0(4) universality class in 3 dimensions. The singularity which deter- 
mines the radius in (^Q) lies in the complex plane and pinches the real axis (together with its 
conjugate) at the critical point. The convergence radius /x_r(T) has a certain nonanalyticity 
at T = Te, which we shall describe. 

The study of thermodynamic singularities, or partition function zeros, in the complex 
plane was pioneered by Yang and Lee Their analysis was extended by Fisher from 

the complex magnetic field singularities to complex temperature singularities [ii. The 



properties of these singularities following from scaling and universality have been further 
studied by Fisher [3 ; Itzykson, Pearson and Zuber ly| and others. 



From the point of view of the lattice studies of the QCD phase diagram, we would like to 
distinguish two separate issues. One is the convergence of the series as the truncation order 
is increased. This is the issue which this study will impact. The other is the convergence of 
each term in the series to its thermodynamic limit as the volume and/or the number of Monte 
Carlo configurations are increased. The sign-problem is affecting this latter convergence and 
will not be addressed here (see, e.g., 171 llSl. llOj). 
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FIG. 1: A sketch of the QCD phase diagram in the vicinity of the critical hne HciT) for zero and 
nonzero quark mass. 



II. UNIVERSAL PROPERTIES OF THE COMPLEX SINGULARITIES 

Here we describe generic universal properties of the complex thermodynamic singularities. 
The basic results following from scaling and universality are not new, and can be found in 
[l6l |. For clarity and completeness, we rederive the needed facts here using slightly different 
approach. Our purpose is to apply these results to QCD at finite temperature and density. 



A. Complex /j, singularities as Fisher zeros 

We shall first consider QCD in the chiral limit - the limit of two lightest quark masses 
taken to zero. The phase diagram in the (T, fi) plane is sketched in Figure H The low- 
temperature phase is separated from the high-temperature phase by a phase transition. This 
is unavoidable, because the symmetry of the ground state must change from SU(2)yxU('l)B 
to the full symmetry of the action SU(2)y xSU(2)yixU(l)B as the temperature is raised ^20i] . 
This transition is of second order for fi < fi^, and of first order for fi > fi^. 

Let us now fix T at a value in the interval (T3, Tc) and study the behavior of singularities 
in the complex fi plane. On the real axis, as fi increases from zero, the transition occurs at a 
value which we denote as ^c{T) (see Fig. P). Therefore, at a fixed temperature, the change 
of the chemical potential /i — fic{T) is a relevant perturbation. In the universality class of 
the 0(3)— 5>0(4) transitions, to which QCD chiral restoration transition belongs j2l|, there 
is only one relevant variable - the thermal variable t (we are discussing the symmetry limit 
- the magnetic field variable h is absent). Therefore, we are led to consider the universal 
behavior of the singularities in the complex temperature plane, which correspond to Fisher 
zeros. The scaling parameter t is to linear order 

t~/i2-/i,2(r). (2) 

We use /i^ instead of /i to ensure that at /i = the thermal variable is proportional to T — T^. 
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B. Partition function zeros and electrostatic analogy 



As first exposed by Lee and Yang [13[, the thermodynamic singularities in the complex 
plane are related to the zeros of the partition function. For a finite system the partition 
function Z, by definition, is strictly positive for real values of the parameters. However, Z 
has zeros in the complex plane, whose number grows linearly with the size of the system. In 
the thermodynamic limit the zeros typically coalesce into cuts. A phase transition occurs 
where such a cut pinches (second order) or crosses (first order) the real axis. 

The considerations of Lee and Yang apply to the variable A = exp(/i/T) since the partition 
function of QCD is a polynomial in this variable due to quantization of the baryon charge. 

It is convenient to use the electrostatic analogy. The partition function - a polynomial, 
can be written in terms of its roots 

Z{\) = n(A - A.) . (3) 

k 

Therefore the free energy (or grand potential, if we are dealing with grand canonical ensem- 
ble, as in QCD) is 

^](A) = -TlogZ = -T5^1og(A - \k). (4) 

k 

For complex A, the real part Refi can be interpreted as an electrostatic potential created 
by charges located on the plane (Re A, Im A): 

Re(](A) = -T^log|A- Afc|. (5) 

i 

All charges have the same magnitude and sign (degenerate roots can be treated as coincident 
charges). Now let us assume, as is known to be true in most cases, and in the universality 
region of interest in particular, that the zeros coalesce into 1-dimensional curves in the 
thermodynamic limit. Then, the electrostatic potential Refi is continuous across such a 
curve, while the analog of the electric field 

£;=-V(ReSJ) = -(— .— j=(-Re.Im)- (6) 

is discontinuous - the normal component jumps by an amount proportional to the linear 
charge density p on the curve. This curve can be viewed as the location of a cut on a 
Riemann sheet of the analytic function Q{X). 

C. Stokes boundaries in the scaling region at /i = 

Consider now the critical region in the vicinity of a critical point Ac. In the scaling regime 
the singular (i.e., non-analytic) part of the potential Q{t) is proportional to a power of t: 

n- t = (7) 

i^smgW - I (_^)2-a^ t<0; Ac 

which defines the universal specific heat exponent a and the amplitudes A±, whose ratio is 
also universal. Off the real axis fl{t) must be an analytic function everywhere except for 
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FIG. 2: Universal behavior of Stokes boundaries in the scahng region for zero and non-zero sym- 
metry breaking parameter h. Only upper complex half-plane is shown. The trajectory of the 
branching point t^{h) is indicated by a dashed line. 



discontinuities across the cuts. Such cuts (at least two, by symmetry Q{t*) = ^l{t)) must be 
present, because the function ((7j) taken at t > does not match its analytic continuation from 
the t < axis along a path around t = 0. The location of the cuts can be determined using 
electrostatic analogy, which requires Kefl to be continuous across the cut. Parameterizing 
t = —s e*"^ using real parameter s > we find: 

A+cos[(2 - - vr)] = A_ cos[(2 - (8) 

Therefore, the cuts are straight lines at an angle with respect to the negative t axis given 
by (cf. 0) 

, , cos(7ra) — y4_M+ , . 

tan 2 - a)ip] = ^ / / + , 9 

sm(7ra) 

as shown in Figure |21 All quantities entering this formula are universal. 

The cuts are termed Stokes boundaries in jl6| — they carry conceptual resemblance to 
the anti-Stokes lines in the WKB theory. Across these lines, the function Q switches from one 
of its Riemann sheets to another. The density of the "charges" on the cut is proportional to 
the discontinuity of the normal component of E, and thus to Im (e**^ dfl/dt), which vanishes 
at the branching point as 

p~|t|^-". (10) 



D. Stokes boundaries at /i / 

The magnetic field h is another relevant variable near the 0(4) critical point. The mag- 
netic field breaks the 0(4) down to 0(3), and in QCD this role is played by the quark mass 
m (more precisely, the average of the u and d masses). At h ^ the free energy is analytic 
function of t at t = 0. In the scaling region the singular part of the free energy scales as a 
power of h if t is also changed to keep the scaling variable x = th~^^^^^^ fixed. Therefore 
the (two) branching points must be located at a point away from the origin given by 

U = x^h"'^"'^ (11) 

and at (t*)*, where x^, is a complex constant. The phase of t,, (the polar angle coordinate 
of the branching point) is determined by the following argument. By scaling postulate, the 
singular contribution to the free energy must be given by 

aing = /^'+'/'A(t/l-l/(^^)) (12) 
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where the function A{x) is analytic at a; = 0. A{x) has two cuts which originate at points 
X* and (x*)* and go off to infinity. Consider fising as a function of complex h at fixed 
real t. By symmetry {h ^ —h, h ^ h*), the Stokes boundaries of this function lie on the 
imaginary axis (more rigorously, this follows from the Lee- Yang theorem and universality). 
Thus the branching point, for t > 0, at /i* = {t/x^Y^ , is purely imaginary, and therefore 
PSaigx^ = Tx/2. Thus, 

U = \x^e^^h^'^^'\ ^ = -^. (13) 

To summarize, at /i = 0, the complex singularities in the scaling region of the thermal 
parameter t form cuts (Stokes boundaries) which go along the rays at angle ip with the neg- 
ative real axis given by With increasing symmetry breaking parameter h the branching 
point shifts away from t = by an amount proportional to h^/^f^^'> along the direction at the 
angle to the positive t axis given by (fT^ . This is illustrated in Fig. |21 



III. SINGULARITIES OF QCD IN THE COMPLEX PLANE 
A. Chiral limit: m = 

Let us begin with the chiral limit m = 0. Consider the interval of T G (T3,Tc). In this 
interval increasing /i leads to the second order transition at /z = ^c(T). As discussed in the 
previous section in the vicinity of the transition we must identify t with fi"^ — fic^{T) (up to 
an irrelevant constant factor) - Eq. Q. Thus, at a given temperature T, near ndT) the 
location of the singularities in the complex /i plane is determined by the universal arguments 
of the previous section. Conformal transformation /i — ^ does not affect the angles and 
ip away from /x = 0. 

In particular, at m = 0, the (two) cuts should originate at the branching point located 
at fJ,c(T) on the real axis, and follow the rays at angle ip given by Q with respect to the 
negative real axis as shown in Figure |21 (left). 

Taking the values a ~ —0.25 ^2] and A^/A_ ^ 1.6 [23] we estimate the value of the 
angle as ~ 77°. At the tricritical point a = 1/2 and A^/A_ = and thus ip = 60°. 



B. Small quark mass: m ^ 

At finite quark mass m the second order line /i = Hc{T) is replaced by an analytic 
crossover for all temperatures T > Te- The critical ending point of the first order transition 
is located at a temperature we denote Te- For T < Te the transition is of the first order. 

At fixed T > T^, and small mass m, in the vicinity of the crossover point the singularities 
are described by the universal arguments with 

/i~m. (14) 

For the purpose of the discussion we can define the crossover point as the value of the 
real part of the branching point: /^crossover = Re (/i*(m)). The branching point, /i*(m) or the 
nearest singularity to the real axis, is shifted by the amount proportional to m^^^f^^\ which 
is m°'^^, using the 0(4) critical exponents p^], in the direction along the ray at angle (fT^ 
Ip = 7r/(2/55) ^ 48° towards the positive real axis as shown in Figure IHl (left). 
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FIG. 3: Stokes boundaries in QCD at fixed T and two different values of m (left) and at fixed 
small m and three different values of T (right) as dictated by universality. 



When T is decreased towards Te this branching point ^^{m) (and its conjugate) approach 
the real axis and pinch it when T = Te- At this temperature, the cuts originate from the 
branching point on the real axis ^ie (see Figure El (right)). The point {Te.^ie) is the QCD 
critical ending point. This ordinary critical point is in the universality class of the Ising 
model 0, 1^ . The initial direction of the cuts near this point is perpendicular to the real 
axis, i.e., ip = 90°. This follows from the fact that the perturbation fi — fiE is magnetic-field- 
like, 

hr-^ H- He, (15) 

and from the fact that singularities in the h plane lie on the imaginary axis. 

The reason that /i — fiE is not t-like, but /i-like, is the following. In the vicinity of the 
critical point, since the 0(4) is explicitly broken, the perturbation h — he affects both the 
thermal variable t as well as magnetic- field-like variable h to linear order. Since (35 > 1, 
the scaling variable x = th~^^^^^^ is small, which means the perturbation fi — he takes the 
system into the region where the variable h dominates the scaling.^ 



IV. RANDOM MATRIX MODEL 



In this section we illustrate the universal properties discussed above by determining the 
complex plane singularities of a random matrix model of the QCD partition function at 
finite T and /x which was introduced in Ref. and applied to the study of the QCD phase 
diagram and the (tri) critical point. At yU = this model is equivalent to the finite-T model 
of Ref. 2^1 and at T = - to the finite-/^ model of Ref. ji^. The parameters T, /i and m 
used in this section are dimensionless, and correspond to measuring T in units of = 160 
MeV, fi in units of 2.27 GeV, and m in units of 100 MeV. For more details, see Ref. 



An equivalent way of saying this is by comparing the scaling dimensions of the variables coupled to thermal 
and magnetic relevant operators: yt — l/i' and yh — PS/v. Since both operators couple linearly to the 
variable /i — /i^, and yh > yt, the magnetic field operator dominates the response to jji — perturbation 
near the critical point. In contrast, at m = 0, the coupling to magnetic operator is forbidden by symmetry. 
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The model can be solved in the thermodynamic limit, which corresponds to the infinite 
size of the random matrix, oo, by using the replica trick. This gives, for real /i, T 

and m, 

log Zrm = - min (16) 

<p 

where 

= 02 - i In { [(0 + - (/i + zT)2] . [(0 + - (/. - ^T)'] } . (17) 

Analytically continuing into the complex /i-plane, one finds the branching points of the 
partition function by solving a system of two algebraic equations 

TT = U; = ( branchmg pomts) (18) 

d(p \d<p) 

for two unknowns: and /x. The second equation states that two of the solutions determined 
by the first equation are coalescing into one at this value of /i. The Stokes boundaries can 
be determined by solving the condition Ref2(0i) = Ref2(02) where 0i and 02 are the two 
solutions of the first equation in (|18p which coalesce at the branching point: 

0; Re n(0i) = Re fi(02) . (Stokes boundaries) (19) 

At finite N the partition function can be written explicitly as a polynomial (apart from 
an irrelevant constant factor) 



N/2 



4m= E {^^{^yN-k^-k2)\iF,{k^ + k2-N-l--m'N) 



(20) 



using the procedure similar to [26|, where i-Fi(a; 6; c) is the Kummer confiuent hypergeomet- 
ric function (in ()2U|) it is a polynomial in m?). The zeros are found numerically for N = 120 
and plotted in Figure El together with the Stokes boundaries given by (fT^ . 

Near the point T = Tc, fi = 0, the thermal scaling variable t is proportional to (T — 
Tc) + Cfi^, where C is the constant giving the slope the second order transition curve (see 
eq. (121) )• Therefore it is convenient to plot the zeros in the complex /i^ plane. 

Only the vicinity of the origin /i^ = is shown in Figure 01 The universal properties 
described in the previous section are manifest. It is clear from the form of the solution 
that the critical exponents near the second order line have their mean field values, and 
correspondingly, the angles are ip = 45° and ip = 60°. At the tricritical point, the exponents 
are given by their mean field values also in QCD (albeit with logarithmic corrections): 

= 60°, ip = 72°. One can also see that the density of the zeros decreases near the 
branching point, as dictated by (fTUj) . 

V. CONVERGENCE RADIUS 

Relevant for the search of the QCD critical point is the question of the convergence 
radius /^^(T) of the Taylor expansion around = of the QCD thermodynamic potential 
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FIG. 4: Stokes boundaries and zeros of the N = 120 random matrix partition function at repre- 
sentative values of T at zero and nonzero quark mass. The trajectory of each branching point as 
a function of T is indicated by a dashed hne. 
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FIG. 5: The convergence radius clS cl function of T in the random matrix model at m = 0.07 
(7 MeV). The value of is the distance of the singularity on the dashed line on Fig. |l]from the 
origin. The critical point at T = Te where the singularities pinch the real axis is shown. 

as a function of T. For sufficiently small quark masses, and sufficiently near Tc, the position 
of the nearest singularity, limiting this radius, is determined by the universal arguments 
given above. For illustration, this radius in the random matrix model is plotted on Figure |31 
All generic and universal features are manifest in Fig. As T decreases, and the branch- 
ing point singularity slides along the dashed line on Fig. I^from left to right, the radius 
contracts, and then, below the crossover temperature, begins to expand again. From the 
universality arguments of Section ITTTl (see Fig. El (right)) we conclude that near the chiral 
(m 0) limit the minimum value of the radius scales with m as 

mm filiT) ~ mi/(^^) ~ m°-^^ (21) 
and is achieved at a temperature T which scales as T — ~ m^^^^^'' ~ mP'^'^. 
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Further away from the minimum, at the critical point, T = Te, the singularity and its 
conjugate pinch the real axis. At this point one observes a non-analyticity in /i|.(T): a 
contribution of order {Te — T)^^ turns on below T = Te- This is due to the kink in the 
dashed line on Fig. HI at Te- In QCD P6 ~ 1.56, given by the exponents of the 3d Ising 
model, while in the random matrix model (36 has the mean field value 3/2. More explicitly, 
the trajectory of the branching point fi^ near fiE is given by: 

/i,(T) =fiE + Ci{Te -T) + i C2{T - Te)"' + O {{T - Te)') , (22) 

with some nonuniversal positive coefficients Ci_2- To derive this equation one observes that 
both t and h scaling variables are linear combinations of (T — Te) and {n — he) and uses 
equation (|13|) for the branching point. The fact that the third term on the r.h.s. in (|22p 
is purely imaginary for T > Te is related to the fact that the branching point in the h 
plane is on the imaginary axis for t > as discussed in Section IIIDl Therefore 

f^UT) = |/i*(^)r = /^l + ci{Te -T) + C2 e{TE - T){Te - T)"' + O ((T - Te)') - (23) 

Below Te, the singularity continues to move away from the origin, and the radius of 
convergence continues to increase. The radius is now determined by the spinodal point of 
the first order phase transition (but see discussion in Appendix). This singularity resides on 
the continuation of the physical Riemann sheet under the cut.^ 

One must point out that the random matrix model of Ref. Q does not capture a known 
feature of the QCD partition function - the periodicity, or invariance under the shift /i — »• 
fi + 27ciT, which is due to the quantization of the baryon charge.^ As shown by Roberge 
and Weiss [i^ this periodicity is related to the appearance of a Stokes boundary given by 
Im /i = ttT for sufficiently high temperatures. For T of order 160 MeV, this Stokes boundary 
could interfere with convergence of the series only if the singularity we discuss moves further 
than ^ 500 MeV. 



VI. SUMMARY AND DISCUSSION 

We have described the location as well as temperature and quark mass dependence of 
the singularities of the QCD partition function in the complex n plane. In the vicinity of 
the chiral phase transition at m = the universality and scaling arguments predict that in 
the infinite volume the singularities are two complex conjugate branch cuts originating at a 
branching point on the real fi axis. The cuts are oriented at an angle to the negative /i axis 
given by Q. At nonzero m the branching points (and the cuts) are shifted in the direction 
given by angle ip ~ 48°, by a distance of order m*^'^*^ (see Fig. |2)). 

A related consequence of the universal behavior of the complex T singularities, and the 
fact that ip < 90°, is the prediction that the crossover point at m 7^ 0, defined as the 
projection of the closest singularity onto the real axis, is above the second order 0(4) line, 
as sketched in Fig. [TJ 



^ On a more subtle level, one has to note that the fact that the singularity H22|l remains on the real axis 
is an artifact of the mean field critical behavior in the random matrix model {/3S — 3/2). In QCD, the 
singularity moves off the real axis by an amount which scales as {Te —-T)^^. 
A random matrix model which docs capture this feature is studied in |27j . 
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The singularities we describe determine the convergence of the Taylor expansion around 
the point = 0. As a result, the radius of convergence at T = Tc is limited by a 
singularity whose distance from the origin scales as fiR^ ~ m°'^^, vanishing in the chiral 
limit. At T = Te the convergence radius fJ.R(T) shows nonanalyticity described by Eq. ((221) • 
The random matrix model of Ref . 01 illustrates these universal predictions: Figures ^ El 
The knowledge of the complex plane singularities might be used to improve the Taylor 
expansion methods, for example, by constructing Fade or similar extrapolations ac- 
commodating the correct universal singular behavior. It can also be used to crosscheck the 
results of lattice Monte-Carlo simulations, by comparing the expected universal behavior of 
the partition function zeros to the output of a lattice calculation. 
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APPENDIX: TAYLOR EXPANSION NEAR A STOKES BOUNDARY 



The following question has the practical importance for lattice Taylor expansion methods 
3 Does a Stokes boundary limit the convergence radius of a Taylor expansion and 
how? A view often expressed in the literature is that the thermodynamic functions can be 
analytically continued past the Stokes boundary (or a first order phase transition) - the 
true singularities being located at the branching points, where the Stokes boundaries end. 
Although this is correct in the strictly infinite volume V = oo, this certainly is not correct 
for any finite volume, no matter how large. This must be so since the singularities of the 
free energy of the type log(z — z^) are located along the Stokes boundaries. How then do 
these singularities appear not to limit the radius of convergence in the — ^ oo limit? The 
answer is that they do, but in order to see the divergence in the series, one needs to continue 
the expansion beyond some large order which slides to infinity asV^oo. We determine 
n* in this Appendix. 

Consider a simplified problem, where the expansion point zq is located sufficiently close 
to the Stokes boundary so that the curvature of the boundary can be neglected, and the 
zeros can be considered equally spaced and extending to infinity. It is also convenient to 
transform the variable z to conformally map the Stokes boundary to the imaginary axis as 
shown in Figure IHl 

The contribution to the thermodynamic potential from the Stokes boundary at finite 
volume is the given by ^ 



n 



sing 



loghv^^n 



k=l 



Txk 



log 2 siYihiyzV) = -zV . 



(A.l) 



^ A more common example il = — log 2 cosli(zl^) can be analyzed similarly, with a little additional compli- 
cation, unnecessary for this discussion. 
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A 



FIG. 6: The convergence disk of the Taylor series around point zq near a hne of complex plane 
singularities. 



We have also rescaled the variable z so that the spacing between zeros is vr/l^, to emphasize 
the fact the the density of zeros grows linearly with V . In the thermodynamic limit 



V 



-z, 

z + ivr, 



Re z > ; 
Re 2 < . 



(A.2) 



On the Stokes boundary, Re z = 0, Re is continuous but the normal derivative of Re f2 
is not, in accordance with the electrostatic analogy ©. However, the function —z has no 
singularity on the Stokes boundary and can be analytically continued through it. Let us 
choose Zi^ to be on the positive real axis and Taylor expand around such a point (see Fig. P): 



^sing('2^) ^sing(^2^0 



(Zo) = -{z- Zq 



oo ^ 



zq- Z 
Zq 



where 



^ oo 

n! ^-^ 

p=i 



■2zoVp 



(A.3) 



(A.4) 



For large V and n the coefficients approach^ 



n <t:n^ = 2zoV ; 



(A.5) 



We see that the late terms n > n^, ^ V will cause the divergence of the Taylor series outside 
the circle of radius \zq\. The value of n^, can be written in terms of two quantities: the 
distance d of the expansion point from the Stokes boundary and the spacing S between the 
zeros at the given volume (shown in Fig.lH}. Recalling that in terms of the rescaled variable z 
the spacing is 6 = n/V and the distance is d = \zq\, we conclude: 



27rd 



(A.6) 



^ In the interval n.^, <^ n nj/(27r) the coefficients oscillate with maxima reaching n^l\p2nrrfi before 
settling on a„ = \jn. 
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